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Abstract: We evaluate the couplings of pions in the transitions of positive parity 
heavy mesons, (l + ,2 + ) and (0 + ,l + ), to negative parity ones (0~,1~) using a tech- 
nique which is not limited to the soft-pion limit. This is made through a constituent 
quark-meson model (the CQM model) where the amplitudes are obtained by quark 
loops with mesons on the external lines. The results are in good agreement with 
experimental data. 
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1. Introduction 

Positive parity mesons comprising one heavy quark are the object of a wide experi- 
mental search || and this demands a deeper theoretical understanding of processes 
involving these states. A theoretical description within QCD can be attempted at 
present only by numerical simulations or by the use of QCD sum rules. However 
much of what is known about QCD can also be described by using an effective La- 
grangian, encoding the symmetries and dynamical properties of the theory with a 
limited number of parameters. This is possible because the use of expansion pa- 
rameters, such as the energy or the heavy quark mass, makes sense in some limited 
kinematical regions, whereas in principle the effective Lagrangian description would 
require an infinite number of operators, and correspondingly an infinite number of 
free parameters. The effective Lagrangian approach is usually simple enough to allow 
for analytical calculations, thus giving the opportunity to see in detail the interplay 
between symmetries and dynamical information introduced into the model. 

We are interested here in a constituent quark-meson model (the CQM model 
[0, |3|) in the intermediate range of energies between the confinement scale and the 
1 GeV region. In this model both quark and meson effective fields are present. 
In the following we extend the analysis described in [Q, || to the study of strong 
decays D% — > Dn, D{ — > Dn. Here the states D\, D\ are the observed 2460 and 
2420 charmed positive parity mesons. We also generalise the computation of strong 
decays — > Dir, D — ► D*tt made in avoiding the soft-pion-approximation used 
there (the Df, D states have J p = 1+, 0+). 

Our strategy is to consider a quark-meson Lagrangian where transition ampli- 
tudes are represented by diagrams with heavy mesons attached to loops containing 
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heavy and light constituent quarks. The effective meson-quark vertices emerge as a 
result of path integral bosonization of an underlying 4-quark Nambu-Jona-Lasinio 
(NJL) interaction. Our results are encouraging especially in consideration of the tech- 
nical simplicity that the method offers in comparison to other available approaches, 
such as for instance in |3J] . 



2. The model 

A short presentation of the CQM model is in order; for a full account see || and 
references therein. The model is an effective field theory containing a quark-meson 
Lagrangian incorporating heavy-quark symmetries ||. For our purposes this La- 
grangian can be divided into two parts: 

£ — + C-m- (2-1) 

The first term involves only the light quark degrees of freedom, x, and the pseu- 
doscalar light mass octet. At the lowest order in a derivative expansion one has: 

C u = x(^ M 7„ + A^ l5 ) X - mxx + f -fd^d^ . (2.2) 

o 

Apart from the mass term, Cu is chiral invariant. Here = <9 M — iV^, £ = 
exp(z7r// 7r ), S = £ 2 , f n = 130 MeV. tt is the 3x3 matrix representing the flavour 
SU (3) octet of pseudoscalar mesons and the explicit expressions for V M and A? are: 

A» = - • (2.3) 

As one can see by expanding the exponential of the pseudoscalar mesons £, V 1 
generates couplings of an even number of pseudoscalar mesons to the XX pair, while 
A^ gives an odd number of 7r fields. 

The effective Lagrangian containing heavy and light quarks together with meson 
fields is: 



= Qviv ■ dQ v - 



X(H + S + if^—)Q v + h.c. 



+ ^Tr[(H + S)(H-S)] + ^-Tr[T^] . (2.4) 

Here Q v is the effective heavy quark field of Heavy Quark Effective Theory (HQET) ||, 
A x is the momentum scale characterising the convergence of the derivative expan- 
sion, usually taken as the chiral symmetry breaking scale A x ~ 1 GeV, and H, S, 
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T fields are the effective meson fields discussed in || which, in the explicit matrix 
representation, have the form: 



H 

S 



2 
2 



[p;r - ^7 5 ] 



r 2 IP 



(2.5) 
(2.6) 

(2.7) 



Here P, P* etc. are annihilation operators normalised in the following way: 



(0|P|Q?((T)) 

<o|f*"|W)> 

(0|P 2 *HQg(2+)) 



•/.VI 



H 



M 7 



(2.8) 
(2.9) 
(2.10) 



H, S and T M describe respectively the (0", 1") doublet, the (0+, 1+) doublet D , D*{ 
and the (l + ,2 + ) doublet D\,D\ that are predicted by HQET. The part in the La- 
grangian containing the T field is not directly read from the underlying NJL contact 
interaction, but it constitutes its expected extension. The coupling constants G3, G4 
and the field renormalization constants Zg, Zt are discussed in || (for numerical 
values see below). Another parameter of the model is Ah, i.e. the difference between 
the H meson doublet mass and the mass of the heavy quark involved. This parameter 
is not predicted by the model but, once Ah is fixed, the model gives As (as a result 
of the dynamical information introduced into the Lagrangian by the underlying NJL 
contact interaction), while is still a free parameter that can be estimated using 
experimental information. The following table was computed in ||: 

We allow for a variation of Ah = 0.4 ± 0.1 GeV, since 
only this range of values is allowed by the results on semilep- 
tonic weak decays (for a discussion see @). 

The last step necessary to the definition of the model 
is the choice of a regularization procedure. We use the 
Schwinger proper time regularization method assuming as 
UV and IR cut-off A ~ 1.25 GeV and /x ~ 0.3 GeV respec- 
tively. The former ensures the suppression of large momen- 
tum flows through light quark lines in the loops, since the heavy quark carries most 
of the energy and momentum in the system; the latter, being of order of Aqcd, drops 
the unknown confinement part of the quark interaction (see for example 0). The 
regularized form of the heavy quark propagator is: 



A H 


A s 


A T 


0.3 


0.545 


0.635 


0.4 


0.590 


0.735 


0.5 


0.641 


0.835 



Table 1: A values (in 
GeV) 



d 4 k f 



k\ + m 2 



d 4 k 



1/A 2 



(2.11) 
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where m is the constituent light quark mass, which we take as m = 0.3 GeV for u 
and d flavours. An alternative way of cutting off high momenta is proposed in ||. 



3. Strong couplings 

As shown in ||, transitions T — ► Hit are described by an interaction Lagrangian 
that in HQET has the the following form: 

C = ^-Ti[HT»{iD^A)ls] + ^Tr[HT^iMh 5 ] + h.c , (3.1) 

X X 

while for the S — > Hit transitions we have: 

C = igTr(HH^ 5 A„) + [ih TrtJIS^AJ + h - c -] ( 3 - 2 ) 
The general amplitudes, for the processes we are interested in, are defined as follows: 

(D + (p')7r-(q)\D*°(p, V )) = i 9l ^q,q u (3.3) 
(D* + (p',e)7r-(q)\D*°(p, V )) = ig*T £V (3-4) 

(D* + (p',e)ir-(q)\D?(p,r))} = ig^^f + g™(q 2 - ^^)] (3.5) 



rriT 



2 2 

m Q — m. 



(D°(p')n + (q)\D+(p)) = m S ^ H (3.6) 
(D (p')7r + (q)\D* + (p,e))=ig 5 q^, (3.7) 

where we are using the HQET definitions of the symbols H, S, T for labelling mesons 
masses. The connection with hi, hi, g, h is: 

h' 

gi = g 2 = 2 v /m H m T —— (3.8) 



/ iTfluTTL'T h' 

* = V — (3 ' 9) 

, h 

g 4 = \fm H m s — (3.10) 

2l7lH /Q 1 1 \ 

95 = —^9 , (3.H) 

In 

where h! = hi + h 2 . 

Our method consists in a loop calculation of a diagram such as in Fig. 1 involving 

a current and two meson states on the external legs. Let us consider the calculation 
of T — > Hit , i.e. the determination of h' . 
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According to the Feynman rules derived from Eq. 
shown in figure may be written as: 



J) and Eq. (O), the loop 



■-yj ZHZTmHrriTq^rf 



, iN c 
16tt 4 



re § d4k Tr[($ + m)7 M 7 5 (/ + m)7 5 (l + $7^] 



(P-m 2 )[(A; + g) ; 



m 2 ](t> • A; + A + ie) 

(3.12) 

The computation of the previous expression can be decomposed in simpler expres- 
sions involving integrals such as as Z, Z a , Z au , and Z uuX (reported in the appendix). 
The numerators of these integrals are the Lorentz structures 1, k a , k a k u , k u k u k x 
respectively. 

Once Eq. (|3.12|) is calculated, the cou- 
pling h! may be obtained by matching the 
previous calculation in the quark-meson 
loop model Eq. ( |2.4|) with the expression 
of the effective Lagrangian, where the only 
degrees of freedom are heavy and light me- 
sons Eq. ( |3 . 1| ) and h! — h\-\-h2- The same 
strategy is followed for the evaluation of 
h and g in Eq. (|3.2|) . For the T — > Hit 
coupling we find: 




Figure 1: Diagram for the loop calcula- 
tion of the T — > Hit coupling h! in the con- 
stituent quark-meson model. 



ti 



ZtZh{ 



2 

m r 



q-K 



I 2 + A T Z(A T ) + —(I 3 {A H ) 
2q* 



/ 3 (A T )) 



+ P(R t (A T ),S t (A T ), q7T )} 



(3.13) 



where the polynomial P(Ri, Si,q n ), expressed as a sum of powers of q w multiplied 
by the linear combinations Ri and Si of Jj and Z integrals, is given in the appendix, 
together with the expressions for Ii and Z. The parameter q n is q n = v ■ q where v 
is the velocity of the incoming meson and g M is the pion 4-momentum. Consistently 
with chiral invariance we neglect the pion mass and we write q^ = (q n , 0, 0, q n ) in the 
decaying particle rest frame. Also, consistently with the effective Lagrangian ( |3.1| ), 
containing the minimum number of couplings, in the evaluation of Eq. ( |3.13|) we 
neglected 0(ql) terms. Two powers of q n have to be kept, since they appear in the 
D-wave matrix element ( ^.3[ )-( P3D . 

Varying Ah in the same range as in Table I, we can evaluate the variation 
of the field renormalization constants Zh, Z$, Zt defining the renormalized fields 
H = H/ y/Z H , S = S/ \fZ~s and T = T / a/^t- We only report values of Z H , Z$, Z T 



for A# = 0.4 which are respectively 2.36, 1.14, 1.07 GeV -1 (see for a complete 
table). 

For the strong coupling S — > Hit we find: 



ZsZh 



R 1 (A s )-R 2 (A i 



i? 4 (Ac 



m 2 Z(Ac 



(3.14) 
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In this calculation we have considered terms only up to the order 0(g*), to be 
consistent with the energy expansion for the chiral fields, as explained before; this 
result goes beyond the soft pion limit approximation (SPL) which is often used in 
these calculations (see the discussion below). Numerically we find: 

h! = 0.65lg;f , (3.15) 

where the error is induced by the variation of Ah in the range 0.3 — 0.5 GeV. For h 
we obtain: 

h = -0.76 ±0.13 . (3.16) 

Here we use the mass determination of the S— multiplet made in i.e. 2.165 ± 
0.05 GeV. The error in this determination does not contribute significantly to h. 
Note that recent preliminary experimental results by CLEO |§ give a mass of 2.461 
GeV for the broad 1 + state attributed to the S multiplet. This difference may be 
understood as a consequence of a 0(l/m c ) correction to the mass formula. 

Using the preliminary CLEO results for the 1 + mass in the S multiplet (it affects 
the pion momentum in (3.14)), we get: 



h= -0.56 ±0.11 (3.17) 



which is consistent with the one obtained by QCD sum rules fTOf : 



h = -0.52 ±0.17. (3.18) 

The approach used in this paper to calculate the strong coupling constants outside 
the SPL is an improvement with respect to 0, where we found h = —0.85 ± 0.02 
using q n — > from the very beginning. For g, the soft pion limit is instead a very 
good approximation since the decay D* — > Dn has a very limited phase space. This 
is confirmed by the evaluation of g made in : 

# = 0.46 ±0.04. (3.19) 



This value agrees with QCD sum rules calculations |TT] and with a result from 



a relativistic constituent quark model ||12|| , giving g = 0.44 ± 0.16 and g = 0.40 
respectively. 

We conclude by calculating the widths of the T and S states to H plus pion, 
where the heavy quark is a c quark, since these states are already observed. These 
two body decay processes, from the chiral heavy meson Lagrangian, were already 
discussed in ||, to which we refer for more details. Using updated experimental 



data for the masses in the phase space integrals, and our determination of the 
coupling h', we obtain: 

h' 2 

T(D*° -»■ D+n-) = 4.59 x 10 7 — MeV = 19.4MeU (3.20) 
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h' 2 

T(D* 2 ° -> D* + ir~) = 1.33 x 10 7 —MeV = 5.6MeV (3.21) 

x 
h' 2 

T(Df -> ZT+tT) = 1.47 x 10 7 — MeF = 6.2MeF (3.22) 

A x 

in which the central value of h' = 0.65 is used and = 130 MeV, A x = 1000 MeV. 
The full one pion widths are a factor 1.5 larger than those in the previous formula 
because of the neutral pion decay channel. As the total D^ width is dominated by the 
one pion mode in the chiral heavy meson Lagrangian, one can use the experimental 
result of 23 ± 5 MeV to extract an experimental value for h! from the previous 
formula. One obtains a value of 0.51, which is in good agreement with the value of 



h' predicted by our model. From ( 3.22 ) one can also deduce the total one pion width 
of .D* , T = 9.3 MeV, which is only half the measured total with 18.9i 3 ; 5 . However 
this state can mix with the broader 1 + state in the S multiplet, the D[ or can be 



affected by significant 1 / m c corrections [ TA \ . Neglecting this mixing we can compute 



the width of the D[ state by using the preliminary CLEO result for the mass §| and 
the computed strong coupling constant h given in fl3.17|) . We obtain 



T(D[ -> D*n) = 240 MeV , (3.23) 

which accounts probably for most of the total D[ width due to the limited available 
phase space. This result can be compared with the preliminary CLEO data for the 
total width, 290t$ 1 ± 26 ± 36 MeV. 

4. Conclusions 

In this paper we have investigated the phenomenological implications for excited 
positive parity heavy mesons of the model proposed in ||. The work is experimentally 
relevant at present for the charm states, for which these transitions are seen, and a 
comparison with theoretical calculations is indeed needed. We have improved the 
technique for diagrams involving a pion vertex on the light quark line, achieving a 
determination of the strong couplings of the pion in the processes S — > Hit which 



turn out to be very close to the QCD sum rule determinations fj"0 |. Moreover, using 
a technique which goes beyond the soft-pion-limit approximation, we have calculated 
the strong coupling of the pion in the transition T Hit, where such a limit could 
not be reliably applied. This coupling is then used to predict, for charmed heavy 
mesons, the decay widths of the T doublet of states 2 + and 1 + to the H states 1~ 
and 0~ plus one pion. 
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A. Appendix 

In the paper we use several integrals that we list in this appendix. 



h 
h 

h(A) 



iN c r TC s d A k N c m 2 m 2 m 2 

1 (k 2 - m 2 ) = (_ ' A 2 "' U> { * 



16tt 4 J (k 2 - m 2 ) 16tt 2 v ' A 2 V 
iN c /" rc g d A k N r „,„ m 2 m 2 N 



J (k 2 -m 2 ) 2 = loVr 2 A 2 ' "u 2 { ' 



16tt 4 J (k 2 - m 2 ) 2 16tt 2 v ' A 2 > 2 

iN c r c s d 4 k 



16tt 4 7 (k 2 -m 2 )(v k + A + ie) 

iN c /■«* d 4 k 
^ >~m^J (k 2 -m 2 ) 2 (vk + A + ie) 

where T is the generalised incomplete gamma function and erf is the error function. 
We also define: 



Z ( A) = i^- f* ™ 

1 } 16W (k 2 - m 2 )[(k + q) 2 - m 2 }(v k + A + ie) 

= C ^ e ' sm2 £ dxeSA2(x) 11 + erf (A( ^ )] (A - 5) 



Where = (q n , 0, 0, q n ) is the pion 4-momentum and A(x) = A — xq w . We can 
observe that the soft pion limit of the preceding expression, i.e. the limit q n — > 0, 
brings Z(A) — > /4(A) as it should be, (cfr. I). We need another auxiliary integral 
and some expressions which are linear combination of the integrals listed. They are 
the following: 

R^At) = m 2 Z(A T ) - I 3 (A H ) (A.7) 

, 2 ry I A \ , (Itv 



R 2 (A T ) = A T Z(A T ) + A? J h (A.8) 

R 3 (A T ) = |(A T / 3 (A T ) - A H I 3 (A H )) (A.9) 
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^4 (A T ) 

Si (A t ) 
S 2 (At) 

S 3 (At) 
S 4 (At) 
S 5 (A t ) 



A 



T 



(J 3 (A T ) - h(A H )) 



e 



J 2 - A T / 2 - A T Z(A T ) 



Ji + A T / 3 (Atf) - m 2 / 2 - m 2 A r Z(A T ) 

, 2 



g^/i + A h / 3 (Ah)) + -y (4(A#) - 7 3 (A T )) 
|/i + ^(/ 3 (A H ) - I 3 (At)) 



^ttA^ 



(A^J 3 (A^) - J 3 (A T )) 



where q n = v ■ q. 

We next define the following polynomial: 



(A.fO) 

(A.H) 
(A.12) 

(A.13) 
(A. 14) 
(A.15) 



P(R h Si, q„) = --— 8^(llmi2i + 4^ - 6S 2 ) + 2g 2 (-176mi? 4 + 14S X + S 2 
HHq w l 

+ 8S 3 + 48S 4 ) + 3^(88mi? 3 + S 3 - 16S 4 - 325 5 ) + 15S 5 ] (A. 16) 
which is used in the equation for hi. 
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